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Abstract — Spike-Timing Dependent Plasticity (STDP) is be- 
lieved to play an important role in learning and the formation 
of computational function in the brain. The classical model of 
STDP which considers the timing between pairs of pre-synaptic 
and post-synaptic spikes (p-STDP) is incapable of reproducing 
synaptic weight changes similar to those seen in biological 
experiments which investigate the effect of either higher order 
spike trains (e.g. triplet and quadruplet of spikes) |1|-|3|, or, 
simultaneous effect of the rate and timing of spike pairs |4| on 
synaptic plasticity. In this paper, we firstly investigate synaptic 
weight changes using a p-STDP circuit |5 1 and show how it fails 
to reproduce the mentioned complex biological experiments. We 
then present a new STDP VLSI circuit which acts based on the 
timing among triplets of spikes (t-STDP) that is able to reproduce 
all the mentioned experimental results. We believe that our 
new STDP VLSI circuit improves upon previous circuits, whose 
learning capacity exceeds current designs due to its capability of 
mimicking the outcomes of biological experiments more closely; 
thus plays a significant role in future VLSI implementation of 
neuromorphic systems. 

I. Introduction 

Numerous studies have shown that timing and patterns 
of spikes play an important role in changing the strength 
of synapses, a phenomenon called Spike Timing-Dependent 
Plasticity (STDP). Presented results in (6) suggest that this 
rule plays an important role in learning and memory in 
the brain. Previous experiments have demonstrated that the 
classical form of STDP which just considers the timing dif- 
ference between pairs of pre-synaptic and postsynaptic spikes 
(p-STDP) fails to reproduce several biological experiments 
including those that consider higher order spike trains (e.g. 
triplets and quadruplets of spikes) |l| |3|. and those that in 
addition to timing difference between pairs of spikes, bring 
the rate of spike-pairs into action of changing the synaptic 
weight [4 |. In 2006, Pfister and Gerstner reported a new STDP 
rule, based upon the triplets of spikes (t-STDP), which better 
replicates the above mentioned experimental outcomes JT]. 
Since timing-based synaptic changes can play an important 
role underlying adaptive changes to the internal connectivity 
structure of spiking neural networks, many studies have been 
conducting physical implementation of these rules in recent 
years @, Q-(TTJ. 

This paper proposes a new VLSI implementation for t- 
STDP This implementation is based on an improved version 
of a p-STDP circuit already presented in J5J which is able to 
reproduce the exponential behaviour of the learning window, 



but it fails to generate some important biological experiments 
presented in [2|-[4j. We show that the proposed t-STDP 
circuit succeeds in reproducing all the experimentally observed 
biological effects with reduced normalized mean square error. 

The remainder of this paper is organized as follows. In 
section [II] p-STDP model and its corresponding circuit are 
presented and described. Section III introduces the new t- 



STDP circuit and describes how it is related to the t-STDP 



learning rule. Section IV provides some information about 



experimental setup including, experimental protocols, data sets 
and error function. Simulation results are provided in section 
|V] Section [VI] compares the proposed design to previous 



works, followed by discussion and conclusion in sections VII 
and[Vln] 

II. Pair-based STDP rule and circuit 

The p-STDP synaptic modification rule governs weight 
changes in synapse based on the timing difference between 
pairs of spikes Fig. |TJa). This weight change dynamic can be 
given as: 
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t prc is the time difference between a 



single pair of post- and pre-synaptic spikes, r + and r_ are time 
constants of the learning window, and A + and A~ represent 
the maximal weight changes for potentiation and depression, 
respectively. 

There are several VLSI implementations that try to mimic 
the dynamics of the p-STDP rule |(5), 0, @- The circuits 
presented in [7], [8| cannot reproduce the required exponential 
behaviour seen in Eq. [TJ while the circuit presented in |5J 
(Fig.[TJa)) can. Since the targeted t-STDP rule in this paper has 
an exponential behaviour, the circuit presented in ]5) would 
be a more suitable candidate for implementing the t-STDP 
model. However, in order to make the circuit in |5) compatible 
with the t-STDP rule, some modifications are needed, leading 
to the circuit shown in Fig. |TJb). The modifications are as 
follows: (i) Since the classical p-STDP model (Eq.[T]i is weight 
independent, in the modified circuit, the weight dependence 
part (shown in the dashed box in Fig. [TJa)) is omitted, (ii) 
Potentiation and depression in the modified circuit are repre- 
sented with increased, and reduced amount of charges stored 
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Fig. 1. (a) p-STDP circuit presented in |5J. (b) Modified p-STDP circuit. (c) Pair-based STDP learning window generated by Matlab simulation and also 
using different transistor process corners for the modified p-STDP circuit shown in part (b) (Similar protocols and time constants to [6] were employed). 



on the weight capacitor, respectively, which is in contrast to 
the circuit presented in B). (iii) Also, in order to simplify 
the circuit, preLong and postLong pulses which should be 
generated by an additional circuitry, were replaced with Vpie 
and Impost. These signals represent the input pre- and post- 
synaptic pulses in the modified circuit, (iv) Furthermore, when 
considering the implementation aspects of the circuit presented 
in |9), using bias voltages for times constants control results in 
significant variation in time constants under various 3a process 
corners. So, in order to make the circuit more robust against 
this condition, the bias voltages were represented as the gate- 
to-source voltage of a number of diode connected transistors 
that are biased by current sources (M3 and M14). Fig. [TJc) 
demonstrate that using this approach results in very slight 
changes in time constants for all different process corners 
when compared with MATLAB simulations. Furthermore, 
these current sources (I po t and idep) can later be used to fine 
tune the time constants when needed. This approach suggests 
that the time constants, as well as the amplitude parameters are 
more robust against process mismatch. In order to facilitate the 
simulation of these circuit, the scaling approach used in similar 
VLSI implementations of synaptic plasticity in |9j, p2[ was 
used, which is microseconds instead of milliseconds, i.e a scale 
factor of 1000. However, in all simulation results presented in 
this paper, the results are scaled back to biological time to 
make the comparison easier. 

III. Triplet-based STDP rule and circuit 

Triplet-based STDP (t-STDP), as its name infers, changes 
the synaptic weight based on the timing difference among 
triplets of spikes. According to the triplet rule presented in Q], 
half of the potentiation and depression interactions in the 
triplet rule are similar to what happens in the p-STDP model. 
It means that if a pre-synaptic spike arrives in a specified 
time window after a post-synaptic spike, the synaptic weight 
should decrease, and it should increase if the reverse order 
of spike pairs happens. However, the second half is where 
the non-linearity of the t-STDP model appears. According to 
the triplet rule, when a pre-synaptic spike happens, it not only 
interacts with its previous post-synaptic spike(s), but also with 
its succeeding pre-synaptic spike(s) as well. The same scenario 
goes for the time a post-synaptic spike happens and it has 
effect on the previous pre-synaptic as well as post-synaptic 



spike(s). 

It is worth mentioning that, the t-STDP model is not just 
a simple change in the degree of freedom of the p-STDP 
model, but it tries to overcome some deficiencies of the p- 
STDP model. In (TJ, it is addressed that the t-STDP model, 
removed two main problems of the p-STDP formula. These 
problems and how the t-STDP solves them are as follows: 
(1) As p-STDP considers just pairs of spikes, for any value 
of A + > 0, if a pre-synaptic spike precedes a post-synaptic 
one, it brings about potentiation, while according to [4], at low 
repetition frequencies, there is no potentiation. In the t-STDP 
model, this deficiency can be solved by setting A 2 + to a small 
value or in the case of the minimal rule to zero. So it makes 
the potentiation very small which can be neutralized by a bit 
of depression, or it can be zero Jl). (2) Considering biological 
experiments in (4), for A t > 0, potentiation will increase 
with the increase in frequency. However, this behaviour cannot 
be generated by p-STDP, as when the frequency of pairs 
of spikes increases, it causes the pairs to interact with each 
other, so it causes no significant potentiation. This problem 
can be solved again by correct tuning the t-STDP parameters. 
In this case, should be strong enough to make the 

potentiation wins over depression and so to have depression 
in high frequencies (Tj. 

Fig. [2] presents the proposed circuit implementation of the 
t-STDP model. In the proposed circuit, there are eight param- 
eters that can be tuned by controlling eight bias currents as 
follows: Idcpi, -fpoti, ^dcp2 and J po t2 represent the amplitude 
of synaptic weight changes for post-pre, pre-post, pre-post-pre 
and post-pre-post combinations of spike triplets, respectively. 
Another control parameter for these amplitude values in the 
circuit is the pulse width of the spikes which was kept fixed 
during all experiments in this paper (1 /is). In addition to 
these amplitude parameters, the required time constants in the 
model for post-pre, pre-post, pre-post-pre and post-pre-post 
combinations of spike triplets, can be adjusted using 7 t( ji, itpi. 
7 t d2 and I tp2 respectively. 

The proposed circuit works as follows: upon the arrival of a 
post-synaptic pulse, V A post ( n ), M2, M8 and M22 switched on. 
At this time, idepi can charge the first depression capacitor, 
Cdepi> through M2 to the voltage of Vd cp i- After finishing 
^post(n)> Kiopi starts decaying linearly through M4 and with 
a rate proportional to 7 t di- Now, if a pre-synaptic pulse, 



Vp re ( n ) arrives at M6 in the decaying period of Vd ep i, namely 
when M5 is still active, the weight capacitor, Cw, will be 
discharged through M5-M6 transistors and a depression occurs 
during the presence of a pre-synaptic pulse in the interval 
of affect of a post-synaptic spike (post-pre combination of 
spikes). Additionally, if a pre-synaptic spike arrives at M13, 
soon before the present post-synaptic spike at M8, the weight 
capacitor can be charged through M7-M8 transistors and 
a potentiation happens. This potentiation happens because 
the current post-synaptic spike is in the time of affect of 
a pre-synaptic spike (pre-post combination of spikes). The 
amount of potentiation depends on Vpoti, which itself can 
be tuned by the relevant amplitude parameter J po ti. Also, 
the activation interval of Mil can be modified by changing 
the related time constant parameter J tp i. Furthermore, another 
contribution to potentiation can occur if a previous post- 
synaptic pulse, V^ ost ( n _i), arrives at M27 soon enough before 
the current post-synaptic happens at M8 and also before a 
pre-synaptic pulse happens at M32 (this is the same pulse 
as for M13). In this situation, the weight capacitor can be 
charged again through M7-M8 and by an amount proportional 
to Vpot2 an d Vp t3- This is a triplet interaction in the proposed 
circuit that leads to the required nonlinearity mentioned in the 
triplet learning rule, appears. A similar description holds for 
the situation when a pre-synaptic pulse occurs at M6, Ml 3 
and M21 transistors. But this time one potentiation and two 
depression events can happen if the appropriate situation is 
provided. 

The first two parts of the t-STDP circuit (on the top left and 
the top right) are identical to the p-STDP circuit presented in 
Fig. |TJb). Also, the two bottom parts of this circuit carry out 
the triplet terms interactions. This circuit is in correspondence 
to the full triplet learning rule presented in Jl) which takes 
into account all four possible potentiation and depression 
interactions. However, as it is shown in following sections, 
only some of the terms are really necessary to reproduce the 
expected biological experiments. This is referred to as minimal 
triplet learning rule in [1] which makes the required circuit 
much simpler and smaller. 

IV. Experiments setup 

In order to validate the functionality of the pair- and triplet- 
based circuits in comparison with experimental data observed 
in biological experiments, experimental protocols, data sets, 
and error function identical to those presented in (Tl, are 
adopted. Four required experimental protocols, two different 
sets of data and the error function are explained in the 
following. 

A. Experimental protocols 

1) Pairing protocol: Pair-based STDP protocol has been 
extensively used in electrophysiological experiments and sim- 
ulation studies fTJ. In this protocol, 60 pairs of pre- and post- 
synaptic spikes with a delay of At are repeated with repetition 
frequency of p Hz. In many experiments p — 1 Hz, however, 
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Fig. 2. Proposed full Triplet-based STDP circuit. 

it has been illustrated in [4] that how altering the repetition 
frequency affects the total change in weight of the synapse. 

2) Triplet protocol: There are two types of triplet patterns 
which are used in this paper. Both of them consist of 60 
triplets of spikes which are repeated at a given frequency 
of p = 1 Hz. The first triplet pattern is composed of two 
presynaptic spikes and one post-synaptic spike in a pre-post- 
pre configuration. As a result, there are two delays between the 
first pre and the middle post, Ati = t pos t — t pre i, and between 
the second pre and the middle post At 2 = t post — t pre2 . 
The second triplet pattern is analogous to the first but with 
two postsynaptic spikes, one before and the other one after a 
presynaptic spike (post-pre-post). Here, timing differences are 
defined as Ati = t post i - t pre and At 2 = t post2 - t pre . 

3) Quadruplet protocol: This protocol is composed of 60 
quadruplets of spikes repeated at frequency of p = 1 Hz. The 
quadruplet is composed of either a post-pre pair with a delay of 
Ati = tposti — tprei < precedes a pre-post pair with a delay 
of At 2 = t post 2 — t pre 2 > with a time T > 0, or a pre-post 
pair with a delay of At 2 = t post2 — t pre2 > precedes a post- 
pre pair with a delay of Ati = tposti — t pre i < with a time 
T < 0, where T = (t pre2 + t post2 )/2 - (t pre i + t post i)/2. 



Identical to [2], in all quadruplet experiments in this paper, 
At = —A*i = At 2 . 

B. Data sets 

The simulations were conducted using two types of data 
sets: The first data set originates from experiments on the 
visual cortex HI which investigated how altering the repetition 
frequency of spike pairings affects the overall synaptic weight 
change. This data set is composed of 10 data points (obtained 
from Table 1 of [1|) that represents experimental weight 
change, Aw, for two different Ai's, and as a function of the 
frequency of spike pairs under a pairing protocol in the visual 
cortex. The second experimental data set that was utilized, 
originates from hippocampal cultures experiments from J2| 
which examined pairing, triplet and quadruplet protocols ef- 
fects on synaptic weight change. This data set consists of 13 
data points obtained from Table 2 of QJ. This data set shows 
the experimental weight change, Aw, as a function of the 
relative spike timing At, Ati, Ai 2 and T under pairing, 
triplet and quadruplet protocols in hippocampal cultures. 

C. Error function 

Identical to [ 1 1 that tests its proposed triplet model simu- 
lation results against the experimental data and reports their 
differences as Normalized Mean Square Error (NMSE) for 
each data set, we verified our circuit simulation results under 
same condition. The mentioned NMSE [ 1 1 is calculated using 
the following equation: 



E = 



1 p 



(2) 



where AwJ lp , AwJ ir and Oi are the mean weight change 
obtained from biological experiments, the weight change ob- 
tained from the circuit under consideration, and the standard 
error mean of Aw l cxp for a given data point i, respectively; 
p represents the number of data points in a specified data set 
(can be 10 or 13). 

In order to minimize the resulting NMSE for a circuit, there 
was a need to adjust the parameters and time constants to 
minimize the resulting NMSE. In the following subsections, 
the circuit simulation results and applied bias currents for 
setting the required parameters, in order to have the minimum 
achieved NMSEs are reported. Both mentioned circuit shown 
in figures |TJb) and [2] were simulated using parameters for 
a 0.35 /im SiGe CMOS process. This process was selected 
because of cost-related issues, rather than the high speed of 
this process. All transistors are 0.7 /im wide and 0.35 /im long. 
The capacitor values are: CV=10 pF, and other capacitors are 
equal to 10 fF in both presented circuits. 

The synaptic weight capacitor in the proposed design oc- 
cupies a large portion of the silicon area, as it is the case for 
almost all synaptic circuits. Therefore, one of our concerns is 
to reduce this capacitor value in order to make the area of the 
design, smaller. Some of our simulation results show that, it 
is possible to scale the capacitor value (and hence, its size) 
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Fig. 3. p-STDP circuit simulation results. Please note, there are no data at 
30 Hz. The values of p, T and At are used as described in (Tj. 

TABLE I 
P-STDP CIRCUIT BIAS CURRENTS 



Data set 


/pot 


^dep 


hp 


ha 


First 


150nA 


150nA 


24pA 


18pA 


Second 


41 On A 


190nA 


20pA 


5pA 



to one/tenth, however, the circuits bias currents need to be 
retuned. 

It should be noted that, during all experiments in this paper, 
the nearest spike interaction, which considers the interaction of 
a spike only with its two immediate succeeding and immediate 
preceding nearest neighbours, was used |1|. 

V. Simulation Results 

A. P-STDP circuit simulation results 

Simulation results for the first data set, which reflects the 
weight change under a pairing protocol and as a function of the 
pairing frequency, p, using the presented circuit in Fig. |TJb) 
are demonstrated in Fig. [3|a). This figure shows how the p- 
STDP circuit which acts based on the standard STDP rule, 
fails to reproduce the observed experimental results in visual 
cortex reported in B). The minimal NMSE obtained in this 
situation was E=7.26 which is consistent with the reported 
minimal achieved error using computer simulation of the p- 
STDP rule in Fig. 6A of [1]. The four required bias currents 
for controlling the model parameters are reported in Table [I] 

Simulation results of the second data set also suggest that 
p-STDP rule fails to reproduce experimental results observed 
in hippocampal cultures. Fig. [3jb) shows the circuit simulation 
results along with experimental results for the quadruplet 
protocol, while Fig. [3jc) and Fig. [3jd) represent the triplet 
protocols for pre-post-pre and post-pre-post combinations of 
spike triplets, respectively. The minimal NMSE obtained in 
this situation was E= 10.76, again consistent with the reported 
results in Fig. 6B of (T). The four required bias currents for 
controlling the model parameters are reported in Table [I] 



TABLE II 

First t-STDP circuit bias currents 



Data set 


-Jdcpl 


A pi 




Jpot2 


-ftp 2 


First 


3()0nA 


40nA 


40pA 


1.5uA 


50pA 



TABLE III 

Second t-STDP circuit bias currents 



Data set 


-Zpotl 


^depl 


Api 


hdi 


fpot2 


Itp2 


Second 


160nA 


130nA 


28pA 


20pA 


400nA 


lOpA 



B. T-STDP circuit simulation results 

In order to test the proposed t-STDP circuit under the 
mentioned protocols and using the two data sets, firstly full 
t-STDP circuit was employed. This circuit is shown in Fig. [2] 
and consists of four distinct parts each of them related to one 
of the pair or triplet combination of spikes. However, as stated 
in QJ, during simulations we observed that only some of these 
combinations are really necessary and play significant roles 
in synaptic weight change under different protocols. So, we 
changed the full t-STDP circuit to minimal t-STDP circuit 
in correspondence to two minimal t-STDP rules in (TJ. In 
these minimal rules, the inconsequential parts of the proposed 
full-triplet circuit are removed to have the minimal circuits. 
As it can be extracted from the last line of Table 3 in (JJ, 
the minimal t-STDP rule which is capable of reproducing the 
expected visual cortex weight change experiments (the first 
data set), set pre-post and pre-post-pre spike combination am- 
plitude parameters to zero. And it means that this rule neither 
require the pre-post interactions of spikes, nor the pre-post-pre 
interactions to take part in synaptic weight modification. So 
we don't need these parts also in the corresponding minimal t- 
STDP circuit. This circuit composed of 19 transistors (exclude 
the parts in the dashed and dotted boxes in the circuit presented 
in Fig. [2]) and it can reproduce very similar results to the 
full t-STDP circuit which contains 34 transistors (Fig. 15}. The 
minimum NMSE obtained for the first data set and using this 
first minimal t-STDP circuit was E=0.64 which is near the 
minimum NMSE obtained by means of computer simulations 
of minimal t-STDP model, E=0.34 (error value obtained from 
Table 3 of (TJ). The five required bias currents for controlling 
the model parameters are reported in Table |H| 

Furthermore, the minimum obtained error for the second 
data set using the second minimal t-STDP circuit is E=2.25. 
The achieved results are shown in Fig. |4|b)-(d). The second 
minimal t-STDP circuit is composed of the whole top parts and 
the right bottom part of the full t-STDP circuit presented in 
Fig. [2] (see the last line of Table 4 in (TJ). The obtained NMSE 
using this circuit is slightly better than the NMSE obtained 
using minimal t-STDP model and by means of computer 
simulations (E=2.9 extracted from Table 4 of (Tl). The six 
required bias currents for controlling the model parameters 



are reported in Table III 



VI. Previous works 

There is a previous VLSI implementation [9|, capable of 
reproducing the mentioned biological experiments (except the 
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Fig. 4. t-STDP circuit simulation results. Please note, there are no data at 
30 Hz. The values of p, T and At are used as described in (Tj. 



quadruplet protocol which has not been shown in (9]) similar 
to the proposed triplet circuit. In terms of functionality, these 
implementations are different, since [9] directly implements 
the BCM rule, requiring the voltage of the neuron to take part 
in learning, and eventually leads to compatible changes in the 
neuron architecture. By contrast, our proposed design is based 
upon triplets of spikes from which we extract the required 
triplet, quadruplet, and pairing frequency experiments. Note 
that our proposed circuit is capable also to reproduce the 
effects of the BCM rule (results not shown but to be presented 
elsewhere). In addition, unlike the circuit presented in [9], our 
VLSI circuit acts as a STDP circuit, which can be simply 
used to connect to other sets of (neuromorphic) neurons 
of choice. Furthermore, we show that the proposed circuit 
not only can reproduce the required behaviour seen in the 
mentioned experiments, but also it can be tuned to mimic 
those experimentally observed behaviour with a small error, 
while [ 9 1 just depicts the behaviour and not the required values 
observed in biological experiments. Also, the proposed circuit 
would require smaller silicon real states and lower power 
consumption when compared to the circuit presented in fJj. 



Besides, our previous work 1 10 1, presented a circuit design 
of t-STDP which reproduced the essential behaviour seen 
in the mentioned experiments except for a non-exponential 
learning window. This exponential learning window is an 
essential feature of the original pair-based and triplet-based 
STDP models presented in |JJ, fT3) . So in the present circuit 
design, the circuits closely mimic the mathematical models of 
STDP rules. This has lead to better synaptic weight modifica- 
tion capability in comparison to our initial reported design. 

VII. Discussion 

The proposed design uses a number of transistors operating 
in the subthreshold regime of operation, so it is expected that 
this circuit will be sensitive to process variations. However, 
the simulation results shown in Fig. [JJc) (a 3cr process corner 
variations) depict that the time constants in the STDP learning 
window are almost invariant and there is a slight variation in 
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Fig. 5. P-STDP circuit simulation results for 1000 Monte Carlo runs for 
mismatch analysis. Each curve represents the weight change for a random 
mismatch between current mirror transistors. The inset figure shown above is 
measured biological data and demonstrates the excitatory postsynaptic current 
(EPSC) and the noisy nature of these data (6J. 

the amplitude of the window when considering these corner 
variations. This can be explained by considering the fact that 
in these simulations similar transistors parameters are assumed 
for adjacent transistors. 

Besides, as the circuits consist of current mirrors that are 
used to replicate the required amplitudes or time constants, 
the effect of transistors mismatch on the circuit performance 
must be considered. Here we report some preliminary results 
concerning the transistors mismatch in the proposed circuits 
using two runs of Monte Carlo (MC) analyses. These analyses 
were conducted on both pair- and triplet-based circuits. In 
these analyses, the threshold voltage was varied by using the 
delvto parameter in Hspice. This parameter was set to an 
absolute Gaussian distribution with a nominal value of zero, 
absolute variation of 26.6 mV for a standard deviation of 3. 
These values were obtained from the MC model parameters 
for the mentioned process. 

The first analysis examines the effect of the current mirror 
transistors mismatch on the p-STDP circuit (Fig. |TJb)) when 
generating the pair-based STDP learning window (Fig. [5J. 
This simulation shows that the design is prone to transistor 
mismatch and cannot reproduce the exact behaviour of the p- 
STDP theoretical model. However, at the same time it demon- 
strates that the total behaviour of Long Term Potentiation 
(LTP) and Long Term Depression (LTD) in the window can 
be regenerated. 

The second analysis was done on the t-STDP circuit, using 
the first data set, in the presence of current mirror transistors 
mismatch. Figure [6] shows how the NMSE for this circuit 
varies for 1000 various mismatch conditions created in a 
MC mismatch analysis. Results for the MC simulation in 
this case shows that the NMSE in the presence of transistor 
mismatch can go almost up to 18. These simulations show a 
significant variation in the NMSE when compared with TT 



model parameter simulations. 

Although the circuit shows a significant variation in the 
NMSE comparing to the NMSE in the absence of variation and 
using the typical model, by employing some digital trimming 
techniques p4| , we can ensure the circuit operation with an 
acceptable value of NMSE. Furthermore, another method for 
curing this variation in the NMSE is retuning the bias currents 
in a way that reduce the NMSE to a minimal point even 
in the presence of device variations. The inclusion of these 
techniques with the proposed t-STDP circuit will be covered 
in a future publication. 

VIII. Conclusion 

This paper presents a new VLSI circuit design for STDP 
learning rule based on triplets of spikes. Simulation results 
demonstrate that the pair-based STDP circuit cannot account 
for complex biological experiments i.e., triplet and quadruplet 
experiments, as well as the effect of pairing frequency increase 
on weight change, while the proposed triplet-based STDP 
circuit can closely mimic the outcomes observed in biological 
experiments. Because of these features, the proposed VLSI cir- 
cuit can play a significant role in future VLSI implementation 
of neuromorphic systems. 
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